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The interplay between spin transport and thermoelectricity offers several novel ways of generating, 
manipulating, and detecting nonequilibrium spin in a wide range of materials. Here we formulate a 
phenomenological model in the spirit of the standard model of electrical spin injection to describe the 
electronic mechanism coupling charge, spin, and heat transport and employ the model to analyze 
several different geometries containing ferromagnetic (F) and nonmagnetic (N) regions: F, F/N, 
and F/N/F junctions which are subject to thermal gradients. We present analytical formulas for 
the spin accumulation and spin current profiles in those junctions that are valid for both tunnel 
and transparent (as well as intermediate) contacts. For F/N junctions we calculate the thermal 
spin injection efficiency and the spin accumulation induced nonequilibrium thermopower. We find 
conditions for countering thermal spin effects in the N region with electrical spin injection. This 
compensating effect should be particularly useful for distinguishing electronic from other mechanisms 
of spin injection by thermal gradients. For F/N/F junctions we analyze the differences in the 
nonequilibrium thermopower (and chemical potentials) for parallel and antiparallel orientations 
of the F magnetizations, as evidence and a quantitative measure of the spin accumulation in N. 
Furthermore, we study the Peltier and spin Peltier effects in F/N and F/N/F junctions and present 
analytical formulas for the heat evolution at the interfaces of isothermal junctions. 
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I. INTRODUCTION 

The central theme in spintronics is the generation and 
control of nonequilibrium electron spin in solids. Un- 
til recently the spin generation has been done by optical, 
magnetic, and, most important for device prospects, elec- 
trical means. ^'^ In a typical device spin-polarized elec- 
trons from a ferromagnetic conductor are driven by elec- 
tromagnetic force to a nonmagnetic conductor. There 
the spin accumulates, with the steady state facilitated 
by spin relaxation. (There are also novel ways to gen- 
erate pure spin currents, without accompanying charge 
currents. ^"^°) The concept of electrical spin injection 
was first proposed by Aronov,^^ and experimentally con- 
firmed by Johnson and Silsbee,^^ who also formulated 
the problem from a nonequilibrium thermodynamics and 
drift-diffusion view.^^'-'^^ An equivalent description in 
terms of quasichemical potentials, convenient to treat dis- 
crete (junction) systems was formulated systematically 
by Rashba.^^ This model, which we call the standard 
model of spin injection, is widely used to describe elec- 
trical spin injection into metals and semiconductors^'^'^ 
and can also be extended to ac currents. 

Until recently one particularly interesting possibility of 
generating spin, by spin-heat coupling, has been largely 
neglected. The generation of nonequilibrium spin by heat 
currents and the opposite process of generating heat cur- 
rents by spin accumulation has already been proposed 
by Johnson and Silsbee^"^ based on nonequilibrium ther- 
modynamics concepts (see also Ref. 17). The spin- heat 
coupling is now the central point of spin caloritronics 
(or spin calorics). ^^'^^ Although the theory of thermo- 
electricity has long been known, ^'^'^-'^ only experimental 



improvements over the past few years have made its ap- 
plication in the context of generating and transporting 
spin appear possible. ^^"^^ 
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FIG. 1: (Color online) Schematic illustrations of the Seebeck 
(a) and spin Seebeck (b) effects. Here AT is the temperature 
difference, V the voltage, j the charge current, js the spin 
current, and vertical arrows denote up/down spin projections. 

At the heart of spin caloritronics is the spin Seebeck 
effect (see Fig. 1).^''"^^ The conventional Seebeck effect, 
also called thermopower, describes the generation of an 
electric voltage if a thermal gradient is applied to a con- 
ductor. In analogy, the spin Seebeck effect describes the 
generation of spin accumulation in ferromagnets by ther- 
mal gradients. The effect was originally observed in the 
ferromagnetic conductor NiFe,^'''^'^ where indication of 
spin accumulation over large length scales (millimeters), 
independent of the spin relax:ation scales in the ferromag- 
net, was found. Since it also exists at room temperature, 
the spin Seebeck phenomenon may have some technolog- 
ical applications.'^^ 

However, the spin Seebeck effect is not limited to 
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FIG. 2: (Color online) Schematic illustrations of the Peltier 
(a) and spin Peltier (b) effects, where j and js denote the 
charge and spin currents. The thermal current jq is different 
in each region. Small vertical arrows denote up/down spin 
projections. 



metals. It has also been observed in ferromagnetic 
insulators^^ as well as in the ferromagnetic semiconduc- 
tor (Ga,Mn)As.'^^ This suggests that the spin Seebeck 
effect does not need to be connected with charge flow. In 
(Ga,Mn)As the sample was even cut preventing charge 
redistribution over the whole slab; the spin Seebeck sig- 
nals were unaffected and in both cases, of compact and 
disconnected samples, the Pt stripes pick up the same 
inverse spin Hall signals. '^^'^''^ The evidence points to a 
mechanism of magnon-assisted spin pumping from the 
ferromagnet into the Pt, producing spin currents there. 
A theory for this spin pumping from a ferromagnetic 
insulator was suggested in Ref. 35. It was predicted 
that phonons can play an important role in the spin 
Seebeck effect, leading to its huge enhancement.^^ Re- 
cent measurements of the spin Seebeck effect in multiple 
(Ga,Mn)As samples also suggest that the spin Seebeck 
effect can be driven by phonons.'^'' In order to explain 
the main trends of the observed temperature and spatial 
dependence of the spin Seebeck effect in (Ga,Mn)As, a 
phenomenological model involving phonon-magnon cou- 
pling was introduced. '^^ 

In addition to the Seebeck effect, there is also another 
thermoelectric effect, the Peltier effect, which refers to 
the evolution of heat across an isothermal junction of two 
different materials due to an electric current being passed 
through the junction. ^"'^"'^ Recently a spin caloritronics 
analog to the Peltier effect, termed spin Peltier effect, has 
been predicted and experimentally observed in a permal- 



loy (Ni8oFe2o)(PY)/copper/PY valve stack.38.39 The spin 
Peltier effect describes the heating or cooling at the in- 
terface between a ferromagnetic and normal conductor 
driven by a spin current (see Fig. 2). 

Another fascinating discovery is that of the ther- 
mally driven spin injection from a ferromagnet to a nor- 
mal conductor. In this experiment thermal currents in 
permalloy drive spin accumulation into copper, detected 
in a non-local geometry.^'^ The structures were of sub- 
micron sizes, so it is plausible that the effects are elec- 
tronic in nature, although magnon contributions to such 
thermal spin injection setups could also be sizable. A 
practical model was introduced in Refs. 38,41 to find, 
with a finite elements numerical scheme, the profiles of 
temperature and spin accumulation in the experimental 
devices. Recently, yet another form of thermal spin flow, 
coined Seebeck spin tunneling, has been demonstrated in 
ferromagnet-oxide-silicon tunnel junctions. Here a tem- 
perature difference between the ferromagnet and silicon 
causes a transfer of spin angular momentum across the 
interface between both materials. 

An important goal for both theory and experiment of 
the spin Seebeck phenomena is to decipher the roles of 
the electronic and non-electronic contributions. It is yet 
unclear under which circumstances the electronic contri- 
bution may dominate. It seems likely that when going to 
smaller, submicron structures in which the spin accumu- 
lation will be a bulk effect, the spin phenomena carried 
by electrons will become important. Similarly, in materi- 
als with strong magnon damping, such that magnons are 
in local equilibrium with the given temperature profile, 
electrons may ultimately carry the entire spin Seebeck 
effect. It is thus important to set the benchmarks for the 
electronic contributions in useful device geometries. This 
is what this paper does: we explore the role of the elec- 
tronic contributions in F/N and F/N/F junctions which 
are subjected to thermal gradients and derive useful an- 
alytical formulas for various spin injection efficiencies. 

Our purpose is twofold: First, we use the drift-diffusion 
framework of the standard model of spin injection pre- 
sented in Refs. 1,2,5 and generalize it to include electronic 
heat transport and thereby derive a theory for charge, 
spin, and heat transport in electronic materials. Sec- 
ondly, we apply this theory to describe F/N and F/N/F 
junctions placed in thermal gradients. While the Peltier 
and Seebeck effects in such structures have been inves- 
tigated in Ref. 43, we focus here on the description of 
thermal spin injection and the investigation of the cor- 
responding spin accumulation. We also look at the spin 
injection in the presence of both electric and thermal cur- 
rents, and find the conditions under which the resulting 
spin current in N vanishes. In all junctions studied we 
present, as general as possible, analytical formulas for the 
spin accumulation and spin current profiles, as well as for 
the thermal spin injection efficiency and the nonequilib- 
rium (spin accumulation driven) spin Seebeck coefficient. 
Moreover, we look at several different setups of the Peltier 
and spin Peltier effects and calculate their respective con- 
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tributions to the heating/cooling at the interfaces in F/N 
and F/N/F junctions. 

The manuscript is organized as follows: Following the 
introduction of the formalism and the basic equations in 
Sec. II, the electronic contribution to the spin Seebeck 
effect in a ferromagnetic metal is discussed within the 
framework of this formalism in Sec. Ill, while Sees. IV 
and V are devoted to the discussion of thermal spin injec- 
tion and related thermoelectric effects in F/N and F/N/F 
junctions respectively. A short summary concludes the 
manuscript. 



II. SPIN-POLARIZED TRANSPORT IN THE 
PRESENCE OF THERMAL FLUCTUATIONS: 
CONCEPTS AND DEFINITIONS 

A. Spin-unpolarized transport equations 

As a first step we will restrict ourselves to the descrip- 
tion of transport in an electronic system which consists 
only of electrons of one species, that is, either of spin up 
or spin down electrons (denoted by the subscript A ='f / J, 
throughout this manuscript). The derivation presented 
here is a textbook matter^*^'^^ and is given here to intro- 
duce the terminology needed for the spin-polarized case 
and to match the concepts from the standard spin injec- 
tion model of Ref. 2. 

If this system is in thermodynamic equilibrium, the 
temperature T and the chemical potential r]{T) are uni- 
form throughout the system. Knowing the chemical 
potential,^^ one can calculate the density of the respec- 
tive electron species under consideration from 



n°[j?(T),T] = J degx {s) fo 



r,{T) 



(1) 



where ks denotes the Boltzmann constant, g\ (e) the 
electronic density of states at the energy e, and fo the 
equilibrium Fermi-Dirac distribution function. Similarly, 
the equilibrium energy density is given by 



[ry(T),T]= J de egx {e) fo 



e - r]{T) 



(2) 



The system is not in equilibrium if an electric field 

—'V(p{x) is present in its bulk. In this case the chem- 
ical potential becomes space dependent. This is taken 
into account by replacing ?7(T) with r](T) + efj,\{x), where 
the quasichemical potential iJ,x{x) now contains the space 
dependence.^'* Since we want to incorporate the effects 
of thermal gradients into our formalism, we furthermore 
allow for different local equilibrium temperatures by re- 
placing the constant temperature T by a space depen- 
dent temperature T{x). As a consequence there is an 
additional position dependence of the chemical potential 
due to the temperature, that is, r/(T) has to be replaced 
by r]\T{x)]. Thus, the total chemical potential is given 
hy i][T{x)]+eijLx{x). 



Assuming the local nonequilibrium distribution function 
to be only energy dependent because momentum relax- 
ation happens on length scales much smaller compared 
to the variation of the electric potential ^{x), one obtains 

hie.x) = fo[ '-'^'^^^\:^'^^^-'^^''^ ]. (3) 

Therefore, the nonequilibrium electron and energy den- 
sities read 



n'\{x)= J de gx{e) fx{e,x) 

nx {V [T{x)] + enx{x) + e^{x),T{x)} . 



ex{x) = y de sgx [s) fx{£, x) 

= el {v [T{x)] + enx{x) + eip{x),T{x)) . 



(4) 



(5) 



The electrostatic field gives rise to an electric current. 

This charge current consists of two parts: the drift cur- 
rent, proportional to the electric field E{x) = —V(fi{x) 
and the diffusion current, proportional to the gradient of 

the local electron density. 

Since the proportionality factor of the diffusion current, 
the diffusivity Dx{e), is energy dependent, it is conve- 
nient to treat electrons with different energies separately. 
The spectral diffusion current density reads 



jDx{x,s)de = eDx{s)V [gxis)fx{s,x)]ds, 



(6) 



from which the complete diffusion current can be ob- 
tained by integrating over the entire energy spectrum. 
The total charge current for electrons of spin A is given 

by 



jx{x) = -axV(p{x) + e j de Dxie)gx{e)V fx{e,x) 



(7) 



where ax is the conductivity. Inserting Eq. (3) into 
Eq. (7), using the Einstein relation, and keeping only 
terms linear in the nonequilibrium quantities Hx{x) and 
(p{x), we find 



jxix) = axW I ''P^ + fixix) } - SxaxS/T{x). (8) 



Here the conductivity is given by the Einstein relation 

<yx = e'^ j de Dx{e)gx{e) (^~^^ ~ e^Dx{eF)g\{£F) 

(9) 



and the Seebeck coefhcient by 
Sx = -— [ deDxie)gx{e) 



-CeT{x) 



dfo 
de 

9xi^F) ^ D'xiep) 
gxi^p) Dxiep) 



r,[T{a 



T{x) 



(10) 
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In both cases the integrals arc calculated to the first 
non-vanishing order in the Sommcrfeld expansion. The 
Lorenz number is £ = {n^ /3){kB /e)^ and g'x{£F) and 
D'x^ef) are the derivatives of the density of states and 
the diffusivity with respect to the energy evaluated at the 
Fermi level sp- 
in addition to the charge current, there is a heat cur- 
rent in nonequilibrium. A treatment similar to that of 
the charge current above yields 

J,,x{x) = 5Aa,r(x)V I^^M ^^^^ 

- CaxT{x)WT{x). 

If the charge and heat currents are defined as in Eqs. (8) 
and (11), currents j\{x) > and jq,\{x) > flow parallel 
to the X direction. 

At sharp contacts the chemical potential and the tem- 
perature arc generally not continuous. Thus, instead of 
Eqs. (8) and (11), discrctized versions of these equations 
are used. The charge current at the contact (C) is given 

by 



"SacSacATc 



(12) 



for the complete local electron density of the system. Ex- 
panding the electron density up to the first order in the 
local nonequilibrium quantities, fJ,^{x), /i^(x), and (p{x), 
and using the Sommerfeld expansion subsequently to cal- 
culate the integrals which enter via Eq. (1), we can write 
the electron density as 



n{x) = no + Sn{x) 



(15) 



Here we have introduced the local equilibrium electron 
density, no = n^r][T{x)] ,T{x)} + nl{r][T{x)] ,T{x)}, 
and the local nonequilibrium electron density fluctua- 
tions, 



6n{x) = eg [iJ.{x) + (p{x)] + egsi^s{x). 



(16) 



Additionally, we have introduced the quasichemical po- 
tential, fi = {fif + /U|)/2, the spin accumulation, /Ug = 
{^f — /i4,)/2, as well as the densities of states g = 
g^isp) + giisp) and gs = g^isp) - gi{ep) at the Fermi 
level. We further assume that there is no accumulation 
of charge inside the conductor under bias (p{x). This as- 
sumption of local charge neutrality is valid for metals and 
highly doped semiconductors and requires n{x) = Hq- 
Hence, Eq. (15) yields the condition 
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and the heat current by 

hxc = TSxe^Xc 0A,7e + A/^Ac) - CTExcAT^, (13) 

where Arjc + eAfixc and AT^, denote the drops of the to- 
tal chemical potential and the temperature at the contact 
respectively. The (effective) contact conductance and the 
contact thermopower are given by "Exc and Sxc respec- 
tively, while T is the average temperature of the system. 



B. Spin-polarized transport equations 

We now consider spin-polarized systems, which we 
treat as consisting of two subsystems, one of spin up and 
one of spin down electrons; each subsystem is described 
by the equations from Sec. II A. 

Energy as well as particles can be exchanged between 
the two spin pools (by collisions and spin-flip processes 
respectively). As energy relaxation (tens of femtosec- 
onds) happens usually on much shorter time scales than 
spin relaxation (picoseconds to nanoseconds), we assume 
that a local equilibrium exists at each position x. Conse- 
quently, both subsystems share a common local equilib- 
rium chemical potential r][T{x)] and temperature T{x). 
On the other hand, the local nonequilibrium quasichem- 
ical potentials iJ,x{x) can be different for each spin sub- 
system. 

From Eq. (4) we obtain 



n(x) =n° {rj [T{x)] + e^i^{x) + eip{x),T{x)} 

+ n° {r, [T{x)] + eii^ix) + eip{x), T{x)} 



(14) 



Sn{x) = 0. 

The local spin density, 

s{x) =n^ {t] [T(,x)] + e/i^(x) + eip{x),T{x)} 

- n° {v [T{x)] + efi^ix) + eip{x),T{x)} 



(17) 



(18) 



can be evaluated analogously to the local electron den- 
sity: First, Eq. (18) is expanded in the local nonequi- 
librium quantities up to the flrst order. The resulting 
integrals arc performed employing the Sommerfeld ex- 
pansion up to the flrst non- vanishing order and, as a flnal 
step, the charge neutrality condition, Eq. (17), is used to 
simplify the result. This procedure yields 



s{x) = so{x) + ds{x), 



(19) 



with the local equilibrium spin density, So{x) = 
n^f^[T{x)],T{x)} 
nonequilibrium spin density. 



nl{r][T{x)],T{x)} and the local 



Ss{x) 



(20) 



It is important to note that So{x) is determined by the 
local temperature T{x), as a result of the rapid energy 
relaxation as compared to the spin relaxation. 

The same procedure can be applied to calculate the 
energy density from Eq. (5), 



e{x) =e° [T(x)] + e^i^ix) + e^{x),Tix)} 

+ el {tj [T{x)] + eii^ix) + e^(x), T{x)} . 



(21) 
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which can be spHt in a local equilibrium energy density, 
eoix) = eO{77 [Tix)] ,T{x)} + el{r, [T{x)] ,Tix)}, and lo- 
cal energy density fluctuations 6e{x), that is, 



e{x) = eo(a;) + Se{x). 



(22) 



Calculating de{x) in the same way as 6s{x), we find that 



Se{x) = 0, 



(23) 



consistent with our assumption of fast energy relaxation 
to the local quasiequilibrium. 

Next, we consider the currents flowing through the sys- 
tem. Since our goal is to calculate the quasichcmical and 
spin quasichcmical potentials, as well as the temperature 
proflle, we not only derive transport equations based on 
Eqs. (8) and (11), but also continuity equations for each 
of the currents considered, that is, charge, spin, and heat 
currents. 

The charge current consists of the electric currents car- 
ried by spin up and spin down electrons, 



j{x) =j^{x)+j^{x) 
,^fv[Tix)] 



(24) 



- -{Sa + S,a,)VT{x), 



where the conductivities are given by u = a| + 174, and 
(jg = (7-|- — (Tj., and the Seebeck coefficients by 5 = S^ + S^ 
and Ss = S-\ — S_i. In nonmagnetic materials (Ts = and 
Ss = 0. In our model we consider a steady state, which 
requires 



Vi(x) = 0, 



(25) 



that is, a uniform electric current, j{x) = j. 

The spin current is the difference between the electric 
currents of spin up and spin down electrons. 



js{x) = jfix) - ji{x) 
1 



+ nix) \ +aVHs{x) 



(26) 



{Ssa + SGs)WT{x). 



As we have seen, the spin density s{x) deviates from its 
local equilibrium value sq{x). Unlike charge, spin is not 
conserved and spin relaxation processes lead to a decrease 
of the local nonequilibrium spin to so(a;). Therefore, the 
continuity equation for the spin current is given by 



5s{x) 



(27) 



where Tg is the spin relaxation time. We will not dis- 
tinguish between different spin relaxation mechanisms 

in our model. Instead, we treat Tg as an effective 
spin relaxation time which incorporates all the different 



spin relaxation mechanisms. We stress that spin relax- 
ation processes bring the nonequilibrium spin s{x) to the 
(quasi) equilibrium value So{x), defined locally by T{x). 
Here we deviate from the treatment given in Ref. 30. 
The heat current. 



jq{x) = 3qA{x) + jq^ix) 

T{Sa + Ssas)^{il[T{x)] 



+ IJ.{x) 



+ '^^^^^±^y^,,{x) - CTaVTix) 



(28) 



is the heat carried through the system by the electrons of 
both spin species. Closely related is the energy current, 



ju{x) =jg{x)- 



V [Tix) 



+ f^{x) > j - lla{x)js{x). (29) 



Inserting Eqs. (24), (26), and (28) and using that the di- 
vergence of the charge current vanishes in a steady state, 
that is, Eq. (25), we find 



Vj„(a;) = [Sj + Ssjsix)] - ^ls{x)Vjsix) 



jCaT{x) 1 



4jC 



VT{x) 



j^{x) jfix) 



fT-|- 



(30) 

where = as /a is the conductivity spin polarization. 
The above formula contains Thomson (first term) as well 
as Joule heating (final two terms). Equation (23) can be 
used to formulate the continuity equation for the energy 
current by enforcing the energy conservation. 



Vi„(a;) = 0. 



(31) 



Thus, if j is treated as an external parameter, the 
transport equation for the charge current, Eq. (24), as 
well as the transport and continuity equations for the spin 
and heat currents, Eqs. (26), (27), (28), and (31), form a 
complete set of inhomogeneous differential equations to 
determine the quasichcmical potentials fJ.{x) and fisix), 
the temperature profile T{x), as well as the currents js{x) 
and jq{x). The solution to this set of differential equa- 
tions, that couple charge, spin, and heat transport, will 
be discussed in the next section. 



C. Spin diffusion equation and its general solution 

In the following the general solutions to the equa- 
tions introduced in Sec. II B will be discussed. Insert- 
ing Eq. (26) into the spin current continuity equation, 
Eq. (27), and using Eqs. (20), (24), and (25) generalizes 
the standard''^ '^^ spin diffusion equation. 



V^lisix) = 



fJ-six) 

A? 



1 



V • [SsVT{x)] 



(32) 
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Here we have introduced the spin diffusion length^'^ 

A« = ^/rsga{l - PDI [e2 {g^ - gf)\. (33) 

As we are primarily interested in Unear effects, we ne- 
glect the position-dependence of the spin Sccbeck coeffi- 
cient Sg, which enters via T(x), and arrive at a simplified 
diffusion equation for the spin accumulation, 



A? 



(34) 



where Ss is evaluated at the mean temperature T. In 
order to solve this equation, we need the temperature 
profile which can be determined from Eq. (31). If only 
first order effects are taken into account, Eq. (31) gives 
the differential equation 



2^s(l-P.2) 



A2 (4£ - S'2 - 5*2 - 2SSsPa 



-n,{x), (35) 



deforming the typically linear profile of T(x). The solu- 
tion to the coupled differential Eqs. (34) and (35) reads 



^,s{x) = Aexp 



Bexp ( -j- 



(36) 



with the modified spin diffusion length 

^'-^i 4£ - (5 + 5.P.)2 • 
Integration of Eq. (24) yields the total chemical potential, 

!?PM + ^(^) = J ^ _ p^^^(^) + ^±|£^T(x) + E. 
e a z 

(39) 

The integration constants A, B, C, D, and E have to be 
determined by including the respective boundary condi- 
tions of the system under consideration. 

If S\ ^ VC (see the next section), it is often possible 
to assume a uniform temperature gradient, that is, 



T(x) = Cx + D. 



(40) 



Then Eq. (34) reduces to the standard spin diffusion 
equation and its solution is given by 

Hs{x)=Aexp(^^^+Bexp(^-^^, (41) 

while integration of Eq. (24) yields the total chemical 
potential, 

^ + K.) = + ^±|^c) . - P„.M + E. 

(42) 



As before, A, B, C, D, and E are integration constants to 
be specified by boundary conditions. However, assuming 
a constant temperature gradient in fcrromagnets is not 
consistent with Eq. (31) and therefore this approximation 
cannot be used in situations which depend crucially on 
the heat current profile (sc!c next section). 

The spin and heat currents can be obtained by insert- 
ing the solutions found above into Eqs. (26) and (28). 



D. Contact properties 

To find the specific solution for a system consisting of 
different materials, such as a F/N junction, we have to 
know the behavior of the currents at the interfaces be- 
tween two different materials. The currents at a contact 
can be obtained by applying Eqs. (12) and (13), giving 



1 



— (S'cSc -|- Ssc^sc) ATc, 



(43) 



(44) 



T (I 

jqc = jqtc + jqic = ^ {Sc'^c + Sgc'^sc) ( + A/X( 

T 

+ — [Sgc^c "I" ^C^Sc) ^I^SC ^TY!iq/^T(;, 

(45) 

where ATc is the temperature drop at the contact, and 
A?7c, A/ic, and A/Xgc are the drops of the local equilibrium 
chemical, quasichemical and spin quasichemical poten- 
tials. Moreover, the contact conductances Sc = S-|-c+^4.c 
and = S-j-c — S^c as well as the contact thermopowers 
Sc = Sfc + Sic and Ssc = S^c — Sic have been intro- 
duced. 

Equations (43)- (45) will be used in Sees. IV and V to 
fix the integration constants of the general solutions, 
Eqs. (40)-(42) and Eqs. (36)-(42) found in Sec. II C. 



III. FERROMAGNET PLACED IN A THERMAL 
GRADIENT 

As a first example we consider a ferromagnetic metal F 
of length L {—L/2 < x < L/2) subject to a thermal gra- 
dient under open-circuit conditions, that is, j = 0. The 
gradient is applied by creating a temperature difference 
AT = T2 — Ti between both ends of the metal which are 
held at temperatures Ti and T2 respectively, as shown in 
Fig. 3. 

At the ends of the ferromagnet we impose the bound- 
ary conditions T{-L/2) = Ti, T{L/2) = T2, and set 
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FIG. 3: (Color online) A schematic illustration of a ferromag- 
net metal placed in a thermal gradient which leads to the 
generation of a spin current. 



js(±L/2) — 0. Since we consider only first order effects, 
the Seebeck coefficients are assumed to be constant over 
the length of the ferromagnet and are evaluated at the 
mean temperature T = (Ti + T2)/2. Using the above 
boundary conditions and Eqs. (36)-(39) yields the spin 
accumulation 

AT sinh(x/A,) 4C-iS + SsPaf 
^'^""^ " 2 i cosh(L/2A3) N{L) 



and the spin current 

3s{x) -- 



5, A, AT 
2 R L 



cosh(a;/As) 



cosh(_L/2As 
4/: - 52 - 5-2 - 2SSsPa 

Ml) ' 

where R = K/ [cj{l - P"^)] and 

N{L) ^AC -S^- Si ~ 2SSsPa 

2^ tanh(i/2As 



(46) 



(47) 



Si i-p 



(48) 



L/2Xs 



If a constant temperature gradient is assumed and the 
reduced model given by Eqs. (40)-(42) is used, the spin 
accumulation reads 



( \ -^\ ^ sinh(x/A,) 
t^s[x)- 2 A, ^ cosh(L/2A,)^ 



and the spin current 

js{x) - 



^ A, AT 
2 R L 



cosh(a::/As) 
cosh(T/2As) 



(49) 



(50) 



where R = A^/ [(t(1 — P^)] is the effective resistance of 
the ferromagnet. 



FIG. 4: (Color online) Profiles of the spin accumulation (a), 
the total chemical potential (b), and the spin current (c) for 
NigiFeig at T = 300 K with L = 100 nm and AT = 100 
mK. The solid lines show the results obtained if a constant 
temperature gradient VT — AT/Lf is assumed, while the 
dashed lines (fully overlapping with the solid ones) show the 
results obtained if the temperature profile is determined by 
V> = 0. 



For metals S\ <C \/~C and Eqs. (46) and (47) reduce to 
Eqs. (49) and (50), that is, the assumption of a uniform 
temperature gradient VT = AT/T is justified. Only at 
the boundaries of the sample both temperature profiles 
differ (insignificantly) as there is a small exponential de- 
cay within the spin diffusion length As ~ As if the full 
model is used compared to a perfectly linear tempera- 
ture profile of the reduced model. 

Equations (49) and (50) from the reduced model cor- 
respond to the profiles of the spin accumulation and spin 
current found in Ref. 47, where a Boltzmann equation 
approach has been used to describe thermoelectric spin 
diffusion in a ferromagnetic metal. 

In Fig. 4 the results calculated for a model NigiFcig 
film with realistic parameters'^*^ [As = 5 nm, a = 2.9 x 10^ 
1/rJm, 5*0 = iS^a^ + Siai)/ia^ + ai) = -2.0 x IQ-^ V/K 
with P^ = 0.7 and Pg = {S^ - S:^)/{S^ + S*;) = 3.0] 
at a mean temperature T = 300 K are displayed. The 
length of the sample is T = 100 nm and the temperature 
difference is AT = 100 mK. As can be seen in Fig. 4, the 
agreement between both solutions is very good. 

Figure 4 (b) shows an almost linear drop of the total 
chemical potential between both ends of the ferromag- 
net. Only at the contacts this linear drop is superim- 
posed by an exponential decay. It is also at the contacts 
that nonequilibrium spin accumulates and decays within 
the spin diffusion length [see Figs. 4 (a) and 4 (c)]. Thus, 
only near the contacts there is an electronic contribution 
to the spin voltage and our electronic model does not re- 
produce the linear inverse spin Hall voltage observed in 
this system, which suggests that a mechanism differ- 
ent from electronic spin diffusion is responsible for the 
detected spin Hall voltage.''^ Also, the "entropic" terms 
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FIG. 5: (Color online) A schematic illustration of a F/N junc- 
tion placed in a thermal gradient. 



in the spin accumulation as introduced in Ref. 30, which 
would lead to a uniform decay of the spin accumulation 
across the whole sample, not just at the distances of the 
spin diffusion lengths off of the edges, do not arise in our 
theory. 



IV. F/N JUNCTIONS 

A. F/N junctions placed in thermal gradients 

In this section we investigate an open (j = 0) F/N 
junction under a thermal gradient. The F/N junction 
consists of a ferromagnet and a nonmagnetic conduc- 



tor, denoted by the additional subscripts F and N in 
the quantities defined in the previous sections. The ex- 
tension of the ferromagnet is given by ~Lp < a; < 0, 
whereas the nonmagnetic conductor is described by val- 
ues < a; < Lat. We also assume that the properties of 
the contact region C, located at x = 0, are known. By 
coupling the F and N regions to reservoirs with differ- 
ent temperatures, T2 and T\ respectively, a temperature 
gradient is created across the junction. The model inves- 
tigated in the following is summarized in Fig. 5. 

Like in the previous section, we can assume uniform 
(but for each region different) temperature gradients 
'^Tp and VTjv and use the simplified spin diffusion equa- 
tion, Eq. (34), and the corresponding solutions, Eqs. (40)- 
(42), to describe the total chemical potential, the spin 
accumulation, and the temperature profile in each region 
separately. The integration constants are solved invok- 
ing the following boundary conditions: T{—Lp) = Ti, 
T{Ln) = T2, and js{—Lp) = js{LN) = 0. Furthermore, 
we use Eqs. (43)-(45) and assume, as in the standard spin 
injection model, ^ that the charge, spin, and heat currents 
are continuous at the interface, giving us five additional 
equations for the integration constants. From this set of 
equations the integration constants, including the gradi- 
ents VTp and VT/v, can be obtained. Depending on the 
choice of the direction of the gradient, one finds that spin 
is either injected from the F region into the N region or 
extracted from the N region by a pure spin current, that 
is, a spin current without accompanying charge current. 

In order to measure the efficiency of the thermal spin 
injection [js(0) < 0] and extraction [js(0) > 0] at the in- 
terface, we calculate the thermal spin injection efficiency 
^ = js{^ = 0)/VT7v, which corresponds to a spin thermal 
conductivity. Our model gives 



(Tat tanh{LN/XsN) {tanh [Lp/Xsp) SscRc (l - Pi 



[1 - cosh-i (Lp/X.p)] SsfRf (1 - P^p)} 



Rp tanh {Ln/Xsn) + Rc tanh (Ln/Xsn) tanh (Lp/Xsp) + Rn tanh {Lp/Xgp) 



r 



(51) 



with the effective resistances for the F, N, and contact 
regions, 

Rn Astv/ctat, (52) 

Rp = X,p/[ap{l-P^p)] , (53) 

R, = 1/ [E,(l - P|)] , (54) 

and the contact conductance spin polarization 

Ps = ^sc/^c (55) 

Equation (51) has been derived in the limit of Sxp/n/c ^ 
\/C, in which the temperature gradients are given by 



AT 



apTZ 



FN 



(56) 



vrAr = 



AT 



(57) 



(58) 



a^TZpN 

where 

7? _ Lp 1 Lm 

l-CpN — r — H . 

Op l^c cfn 

If the sample sizes are large, that is, if Lp ^ X^p and 
Ln ^ XsN, as is usually the case (but not in Figs. 6 
and 7 where Lat < Xsn), the situation at the interface is 
not sensitive to the boundary conditions far away from 
the interface and Eq. (51) reduces to 



(Tn_ 
(7n 



SscRc (1 



Pi) 



SspRp (1 



Rp + Rc 



R 



n 



(59) 



{Ssii^P';))R, 
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where {■■■)r denotes an average over the effective resis- 
tances. The above expressions for the spin injection ef- 
ficiency and the gradients, Eqs. (51)-(59), could have 
also been obtained by using Eqs. (36)-(39) to calcu- 
late the profiles and taking the limit S'af/w/c ^ V^- 
Equation (59) is the spin-heat coupling equivalent of 
the well-known formula for the electrical spin injection 
efficiency. 

Using the spin injection efficiency, Eq. (51) [or Eq. (59) 
for large devices] , the profiles of the spin current and ac- 
cumulation in the N region (0 < a; < Lpf) can be written 
compactly as 



sinh[(a; - Ln)/Xsn] 



js{x) = -kVTn 

and 

which reduce to 

js{x) = KVTArexp(-a;/AsAr) 



sinh(Ljv/Asjv) 
cosh [(a: - Ln)/Xsn] 

sinh(iAr/AsAr) 



and 



(60) 



(61) 



(62) 



(63) 



for Lpf ^ XsN- In particular, at the contact the spin 
accumulation in the nonmagnetic material can be calcu- 
lated as 



//,(0+) = -RnkVTn coth (Ln/X.n) . 



(64) 



Equation (51) also makes it clear that whether there is 
spin injection or extraction depends not only on the direc- 
tion of the temperature gradient, but also on the specific 
materials chosen. 

Another quantity of interest is the total drop of the 
chemical potential across the F/N junction, 

A {rj/e + /i) = [r/(T2) - rj{Ti)] /e + (1(1^) - ^(-^f ), 

(65) 

because — in analogy to the calculation of the total resis- 
tance of the F/N junction in the case of the electrical 
spin injection^ — it allows us to define the total Seebeck 
coefficient S of the device, which can be separated into 
an equilibrium and a nonequilibrium contribution: 



A{r]/e + ^l) = SAT = {So + SS) AT. 



(66) 



Here 



^0 = 



(Sp + SspPaF . 



f {Sc + SscP^) 

2TZfn 



S 



N- 



(67) 

denotes the Seebeck coefficient of the F/N junction in the 
absence of spin accumulation, whereas 



5S 



P^F [M.(-iF) - (0^)] + P^ [a*. (0") - (0+)] 
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FIG. 6: (Color online) Profiles of the spin accumulation (a) 
and the total chemical potential (b) for a NisiFeig/Cu junc- 
tion at T = 300 K with Lf = Ljv = 50 nm and AT = -100 
mK. The sohd fines show the results for i?c = 1 x 10~^^ Q,w? , 
the dashed fines for Rc = 1 x 10~^* Q,vc? . 
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FIG. 7: (Color online) Profiles of the spin current (a) and the 
heat current (b) for a NigiFeig/Cu junction at T = 300 K 
with Lf = Ln = nm and AT = -100 mK. The sohd fines 
show the results for i?c = 1 x 10"^" Q,ro? , the dashed lines for 
Rc^lx 10"" Om^ 



is the nonequilibrium contribution to the Seebeck coef- 
ficient due to spin accumulation. If the extensions of 
the F/N junction are much larger than the spin diffu- 
sion lengths, the nonequilibrium Seebeck coefficient can 
be expressed as 



SS = 



2ap 



TZfn 



(69) 



(68) 



For illustration, the profiles of the total chemical po- 
tential and the spin accumulation are displayed in Fig. 6 
for a junction consisting of NigiFcig (see Sec. Ill for 
the corresponding parameters) and Cu {XsN — 350 nm, 
CTAT = 5.88 X 10^ l/flm, Sn = 1-84 x 10"^ V/K) with 
a temperature difference AT = T2 — Ti = —100 mK 
between both ends of the junction and the mean temper- 



10 



ature T = 300 K.^"''"''^^ Figure 7 shows the spin and heat 
currents for the same system. In Figs. 6 and 7 we have 
chosen = 1 x 10'^^ rjm^ and i?^ = 1 x lO^^^^ rJm^, 
as well as Fs = 0.5, Sc = -1.0 x 10"^ V/K, and 
Ssc — 0.5S'c.^^ There is a drop of the total chemical 
potential across the junction [see Fig. 6 (b)]. For the 
chosen parameters spin is injected from the F region into 
the N region, where nonequilibrium spin accumulates at 
the F/N interface and decays within the spin diffusion 
length [see Figs. 6 (a) and 7 (a) where < \sn]- By 
applying the temperature difference AT into the oppo- 
site direction, that is, by choosing Ti <T2, the situation 
reverses and spin would be extracted from the N region. 
Figure 6 (a) also illustrates that the spin accumulation in 
the N region decreases with increasing contact resistance. 
The heat current flows from the hot to the cold end of 
the junction [jq{x) > 0], as can be seen in Fig. (7) (b). 
Furthermore, one can observe that in the F region the 
heat current is not perfectly constant and decreases at 
X = —Lp as well as at the contact, while in the N 
region the heat current remains constant. 

We now discuss two important cases: transparent and 
tunnel contacts in large F/N junctions where Lp ^ XsP 
and Ljv 3> XsN- For transparent contacts Rc <C Rp, Rn 
and the spin injection efficiency reduces to 



<7n 



SspRp (l — P^f 



2 Rp Rpf ^ ^ 

Thermal electronic spin injection from a ferromagnetic 
metal to a semiconductor, that is, the case of Rn ^ 
Rp, would suffer from the same "conductivity/resistance 
mismatch problem" as the usual electrical spin 

injection does. The nonequilibrium Seebeck coefficient 
can then be written as 

SspXspPaF /, , Rn 



SS = - 



2ap {Lp/ap + Lm/(Jn) 



1 



Rp + Rn 



(71) 

In this case k and 5S are restricted only by the indi- 
vidual effective resistances Rp and i?jv of the F and N 
regions. Moreover, the spin accumulation is contin- 
uous at transparent contacts, that is, ^s{^^) = Ms(0 ) 
and Eq. (64) yields the expression found in Ref. 40 for 

^ls{Q)/^Tp.''^ 

Tunnel contacts, on the other hand, have very large 
effective resistances Rc ^ Rp,Rn for which Eqs. (59) 
and (69) reduce to 

Cat 

" = -T 

and 



Ssc (1 P'P 



(72) 



dS ■ 



2<tf 



S,,(l-P|)[P^Fflf-Ps(fl.F + J?.jv)] 

2 



7^ 



FN 



(73) 

The thermal spin injection efficiency for the tunnel junc- 
tion is determined by the spin-polarization properties of 
the contact and the conductivity mismatch issue does not 
arise in this case. A similar result has also been obtained 
recently in Ref. 52. 




FIG. 8: (Color online) A schematic illustration of a F/N junc- 
tion placed in a thermal gradient with a charge current being 
simultaneously driven through the junction. 



B. Interplay between thermal gradients and 
simultaneous charge currents 

Another interesting effect is the interplay between a 
thermal gradient across the F/N junction and a simulta- 
neous charge current (see Fig. 8). To analyze this process, 
we take Eqs. (40)-(42), this time with a finite charge cur- 
rent J, and replace the boundary condition for the spin 
current at x = —Lp by js{—Lp) = P^pj while leaving 
the boundary conditions for the temperature unchanged 
and also taking js{LN) = as before. By choosing the 
charge current j = jcom appropriately, the effects of the 
charge current and the thermal gradient, each by itself 
applicable for injecting spin into the N region or extract- 
ing spin from it, can cancel each other out. As a result 
we find that for Lp ^ XsP a charge current 



Jc 



Rp(l 



PIf) SsF 



Rc (1 - Pi) 5f 



2npN (RpPaP + RcP^) 



AT (74) 



extracts (injects) the spin injected (extracted) through a 
given temperature difference AT with no net spin current 
in the N region. 

This effect is shown in Figs. 9 and 10 for the 
NigiFcig/Cu junction investigated in this section (see 
above). We find that a current density of jcom = 7.6 x 10'' 
A/m^ (jcom = 1-9 X 10^ A/m^) is needed to compen- 
sate a temperature difference of AT = —100 niK if 
Rc = lx 10-16 nm^ (i?c = 1 X ilm^). Figures 9 (a) 

and 10 (a) show that there is no spin accumulation and 
no spin current in the nonmagnetic material under the 
compensating electric current condition. The drop of the 
chemical potential across the F/N junction is shown in 
Fig. 9 (b) and the heat current flowing from the hot to 
the cold end of the junction in Fig. 10 (b). The spin in- 
jection compensation should be useful for experimental 
investigation of the purely electronic contribution to the 
spin Seebeck effect. 
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FIG. 9: (Color online) Profiles of the spin accumulation (a) 
and the total chemical potential (b) for a NigiFeig/Cu junc- 
tion at T = 300 K with Lf = Ln = 50 nm and AT = -100 
mK if an electric current compensates the spin accumula- 
tion due to the thermal gradient. The solid lines show 
the results for 7?^ = 1 x 10~^® Q.va? , the dashed lines for 



iic = 1 X 10" 



nm^ 



FIG. 11: (Color online) A schematic illustration of a F/N 
junction in the electrical spin injection setup, where (a) refers 
to an isothermal junction and (b) to the situation where 
jq{x) = 0. The fact that in (b) the temperature at one end of 
the junction is not given as an external boundary condition, 
but has to be calculated from the model is implied by "?" . 
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FIG. 10: (Color online) Profiles of the spin current (a) and the 
heat current (b) for a NigiFcig/Cu junction at T = 300 K with 
Lf = Ln = 50 nm and AT — —100 mK if an electric current 
compensates the spin accumulation due to the thermal gradi- 



ent. The solid lines show the results for Rc 
the dashed lines for _Rc = 1 x 10^^'' Qrn^. 



1 X 10" 



nm^ 



Moreover, we remarlc that jcom can be used to describe 
the efficiency of thermal spin injection if one investigates 
an open-circuit F/N junction {j = 0) placed in a thermal 
gradient as above. In this case the spin current at the 
interface, js{x — 0), is described by Eqs. (51) or (59) 
respectively. We can then define the ratio between the 
spin current at the interface and the charge current one 
would have to drive through the junction to cancel the 
thermal spin injection, P = js{x — 0)/jcom- For large 
devices and S 



\F/N/c 



<C V'C this ratio can be calculated 



as 



P : 



{P. 



(75) 



which represents the negative spin injection efficiency of 
the electrical spin injection.^ 

C. Peltier effects in F/N junctions 

As mentioned above, the spin Peltier effect describes 
the heating or cooling at the interface between a fer- 
romagnetic and normal conductor driven by a spin 
current. '^^ In the following we study several different se- 
tups in which a spin current passes through the interface 
of an isothermal (or nearly isothermal) F/N junction and 
which therefore give rise to the spin Peltier effect. 

For every setup investigated in this subsection we as- 



sume L 



N/F 



> A; 



N/F- 



The first setup considered is the 



electrical spin injection in a F/N junction: An electric 
current is driven across an isothermal F/N junction, that 
is, VT = [see Fig. 11 (a)]. Since the entire junction 
is kept at constant temperature, the continuity of the 
heat/energy current, Eq. (31), does not apply and it is 
sufficient to solve just Eqs. (24), (26), and (27), that is, 
the formulas obtained for the electrical spin injection can 
be used. The spin current at the interface is given by^ 



JsW = -F, n , n J = {Pcy)R]- 



Rp ~\~ R(, -f R 



■N 



(76) 



For constant temperature profiles the heat current, 
Eq. (28), is not continuous at the interface and reads 

J,{x)^^3 + ^3s{x), x<0 (77) 



and 



JqW = -^3, x>Q. 



(78) 
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Therefore, the total heat produced/dissipated per time 
at the interface is given by 



where 



and 



rr=J?(0")-J?(o+) = r, + r^ 



T{Sf - S'tv) i 



q 



2 



(79) 


-0.5- 
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(80) 
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(81) 



denote the rates of heat production/dissipation due to 
the conventional (charge) Peltier and spin Peltier effects. 

If the temperature is fixed at just one end of the junc- 
tion, a temperature drop arises across the F/N junction 
due to the heat evolution at the interface. In order to 
estimate this temperature drop, we follow the approach 
used in Ref. 39 and investigate the hypothetical situation 
where no heat enters or leaves the F/N junction and no 
heat is generated inside the junction, that is, jq{x) = 
[see Fig. 11 (b)]. For Sx <^ \fC, the profiles of the chemi- 
cal potential, the spin accumulation, and the spin current 
are nearly identical in the cases of an isothermal F/N 
junction and a F/N junction with jq[x) = (see below) 
and Eq. (82) should give a good estimate for the tem- 
perature difference arising across the junction due to the 
heating/cooling at the interface. 

Thus, instead of VT = 0, we apply the condition 
iq(x) — for any x. This situation requires us to 
solve the full system of differential equations given by 
Eqs. (24), (26), (27), (28), and (31). Since this situation 
depends crucially on the heat current [via jq(x) = 0], 
the full solution given by Eqs. (36)- (39) has to be used, 
which — in contrast to the assumption of constant gradi- 
ents in each region — ensures constant heat currents. The 
temperature far away from the interface is fixed at a given 
value for one region [for example, at T\ in the F region 
as shown in Fig. 11 (b)]. At the interface we impose the 
boundary conditions that the charge, spin, and heat cur- 
rents given by Eqs. (43)-(45) have to be continuous. As 



before, we assume that L 



N/F 



N/F' 



in which case the 



situation at the interface is not sensitive to the boundary 
conditions far away from the interface. Thus, we choose 
lim fisix) = as boundary conditions for convenience. 

a;— ^±oo 

The quantity we are interested in is the temperature 
drop across the entire junction, which can be obtained as 

As usual, S\p/jsi/c ^ y/C and the conventional contribu- 
tion to the temperature drop then reads 



ATch = 



{Sf + SsfPctf) Lf Sc 4- SscPt. 
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FIG. 12: (Color online) Profiles of the spin accumulation (a) 
and the spin current (b) for a NigiFeig/Cu junction with 
Lf = Liv = 1 Atm, 7?c = 1 X 10"^'' ^m^ and j = lO" 
A/m^. The solid lines show the results obtained for an 
isothermal junction at T = 300 K, whereas the dashed lines 
show the results obtained for a junction with jq(x) — and 
r(-Li) = 300 K. 



while the contribution due to the spin accumulation in 
the region around the interface can be obtained from 



AT, = 



S.F (1 - 

SsJl-P? 



2C 



2C 



^''V. (0-) 



(84) 



In this limit the spin current at the interface is given by 
the same expression as in Eq. (76) and we find 

AT, ^^^^^il^^i?^ ((P,)^ - P„p)j 



(1 - Pi) 



+ "^^^ ^' [RpPaF - {Rf + Rn) {Pa)B] j- 

(85) 

In Fig. 12 we display the profiles of the spin accumu- 
lation [Fig. 12 (a)] and the spin current [Fig. 12 (b)] 
in NigiFeig/Cu junctions {Lf = Tjv = 1 /^m, and 
i?c = 1 X 10^^^ rtm^) across which a current j = 10^^ 
A/m^ is driven. As can be seen in Fig 12, the agree- 
ment between the solutions of an isothermal junction at 
T = 300 K and those of a junction where jq{x) = and 
T{~Lf) = 300 K is very good, that is, for Sx y/C 
the behavior of the spin accumulation and current is rel- 
atively insensitive in these cases. 

(82) Having studied the spin Peltier effect in situations 
where the spin current is driven by an accompanying 
charge current, we now turn to a different scenario in 
which we are dealing with a pure spin current (j — 0) 
and there consequently is no contribution from the con- 
ventional Peltier effect. First, we study heating/cooling 

(83) effects at the interface of a F/N junction in the Silsbee- 
Johnson spin-charge coupling setup, ^^'^^ that is, we in- 
vestigate the heat generated at the F/N interface while 
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SxF/N/c ^ a//^ and thus the temperature drop across 
the junction is given by 



AT = 



{Ss{l^Pl))i 
2C 



fis{oo) = -- fis{oo), (88) 



where k is the thermal spin injection efficiency of the F/N 
junction defined in Eq. (59). Equation (88) is the thermal 
analog of the Silsbee- Johnson spin-charge coupling. The 
sign of the temperature drop changes when changing the 
spin accumulation /is(oo) from parallel to antiparallel to 
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V. F/N/F JUNCTIONS 



FIG. 13: (Color online) A schematic illustration of a F/N 
junction in the Silsbee-Johnson spin-charge coupling setup, 
where (a) refers to an isothermal junction and (b) to the sit- 
uation where jq{x) = 0. The fact that in (b) the temperature 
at one end of the junction is not given as an external boundary 
condition, but has to be calculated from the model is implied 
by "?". 



keeping the temperature constant across the entire struc- 
ture, VT = [see Fig. 13 (a)]. The inverse process of spin 
injection, the Silsbee-Johnson spin-charge coupling de- 
scribes the generation of an electromotive force across the 
junction due to the presence of nonequilibrium spin in the 
proximity of the ferromagnet for j = 0. This nonequilib- 
rium spin in the N region generates a spin current which 
then drives the spin Peltier effect. For VT — and the 
boundary conditions fig {—oo) = and /is (oo) ^ (mod- 
eling the spin accumulation in the N region) the standard 
model of electrical spin injection yields 



JsiO) 



Ms (oo) 



Rp + Rc + Rn 



(86) 



for the spin current at the interface. ^'^ Equations (77) 
and (78), which apply to any case of VT = 0, show that 
the heat current vanishes in the N region and the rate of 
heat flowing to or away from the interface is given by 



ptot ps 

9 ^ 9 ^ 



TS 



sF 



Ms (oo) 



2 Rp ~\- Rc -f i?jv 



(87) 



Finally, we look at the Silsbee-Johnson spin-charge 
coupling setup, but instead of keeping the junction at a 
constant temperature, we impose the condition (a;) = 
while keeping one end at a fixed temperature and cal- 
culate the temperature drop across the junction [see 
Fig. 13 (b)]. Applying the additional boundary con- 
ditions lim /is(a^) — and requiring the currents to 

x— f ±oo 

be continuous at the interface, we can use Eqs. (82)- 
(84) with j = 0. Thus, ATch = and the tempera- 
ture drop across the junction is entirely due to the spin 
current/accumulation, AT = AT,. We find that the 
spin current at the interface is given by Eq. (86) for 



A. F/N/F junctions placed in thermal gradients 

The procedure which we used in the previous section 
to describe spin injection in a F/N junction can also be 
applied to more complex structures. Here we will dis- 
cuss spin injection in a F/N/F junction consisting of two 
ferromagnets Fi and F2 (denoted by the additional sub- 
scripts 1 and 2) of lengths Li and L2 and a nonmagnetic 
conductor N (denoted by the additional subscript N) of 
length Ln between the ferromagnets. By adjusting the 
orientations of the magnetization in each ferromagnet in- 
dependently, the junction can be either in a parallel (ft) 
or antiparallel (fi) configuration, that is, we restrict our- 
selves to collinear configurations. The interfaces Ci and 
C2 between the ferromagnets and the nonmagnetic ma- 
terial are located at x = and x — L^- In Ref. 43 
the infiuence of electric currents on the temperature pro- 
file in such structures has been investigated if both ends 
of the device were held at the same temperature. Here 
we consider a different situation: We investigate an open 
curcuit geometry {j = 0) in which both ends of the device 
are coupled to different temperature reservoirs. Holding 
the opposite ends of the device at different temperatures, 
T2 and Ti , gives rise to temperature gradients across the 
junction. Figure 14 gives a schematic overview of this 
geometry. 

The chemical potential, the spin accumulation, and 
the spin current are calculated as in the previous sec- 
tion: Assuming uniform temperature gradients VTi, 
VT2, and VTtv, we use the simplified spin diffusion equa- 
tion, Eq. (34), and fix the integration constants by the 
boundary conditions T(— Li) — Ti, T{Ln + L2) = T2, 
and js{—Li) = js{LM + L2) = 0. Each of the contact 
regions Ci and C2 is characterized by Eqs. (43)-(45) and 
we require that the currents are continuous at each inter- 
face. This allows us to obtain the profiles of the chemical 
potential, the spin accumulation, and the spin current. 

As in the case of the F/N junction, spin is either in- 
jected or extracted at the interfaces between the fer- 
romagnets and the nonmagnetic material. We investi- 
gate the spin injection efficiencies, ki = j,j(0)/VTAr and 
^2 = js(iAf)/VT;v, at the contacts Ci and C2. In gen- 
eral, the expressions for ki and K2 are quite unwieldy. 
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FIG. 14: (Color online) A schematic illustration of a F/N/F 
junction placed in a thermal gradient. 



but can be simplified somewhat if we assume the case of 
Li » Asi and L2 » Xs2- 



Rn coth (Ln/Xsn) + Rcj + R 



Do 



R 



N 



(89) 



^ ^^DosinhiLN/XsN)' 



with i,j = 1,2 and i 7^ j, the thermal spin injection 
efficiencies of the individual F/N junctions, 



Ssr.iRni (1 R'^ 



S.,,.R, 1 - P2 



'2 Ri + Rci + Rn 

as defined in Eq. (59), their effective resistances, 

RpN — Ri + Rci + Rn, 

and 

Da =RIi + {Rci + Ri){Rc2 + ^2) 



(90) 



(91) 



coth {Ln/\sn){Ri+ Rc 



Rc2 + R2 ) R 



N- 



(92) 

Comparing the thermal and electrical^ spin injection effi- 
ciencies of the F/N/F junction, we find that the structure 



of Eq. (89) is similar to the structure of the electrical spin 
injection efficiency. Here the temperature gradient in the 
N region reads 



AT 



(Tat 7?. 



FNF 



where 



FNF 



Li 



1 

^cl 



Ln 

(Tat 



1 

Sc2 



L2 
^2 ' 



(93) 



(94) 



For a given temperature gradient Eq. (89) can be used 
to determine whether there is spin injection [js{0) < or 
js{LN) > 0] or extraction [js{0) > or j,j(Ljv) < 0] at 
the interface C;. The profiles of the spin current and the 
spin accumulation in the N region (0 < x < L^r) are 



jsix) K2 sinh (x/Xsn) ~ sinh [{x — Ln) /X, 



N 



N 



sinh (Ln/Xsn) 



(95) 



and 



fJ-s{x) _ K2 cosh {x/Xsn) - Ki cosh [{x " Ln) /X^n] 
Rn^Tn sinh(I/jv/AsAr) 

(96) 

E Ln ^ XsN, Eq. (89) reduces to Eq. (59), that is, the 
spin injection efficiency of a simple F/N junction. 

In analogy to the procedure employed in Sec. IV we 
can calculate the drop of the chemical potential across 
the F/N/F junction, 

A (r;/e + m) - [riiT2) - r^Ti)] /e + ^i{LN + L2) - M(-ii), 

(97) 

and relate this drop to the Seebeck coefficient S of the 
entire device, 



A {r,/e + n) = SAT = (So + 6S) AT, 



(98) 



which we split into the equilibrium contribution Sq and 
a nonequilibrium contribution SS due to spin accumu- 
lation. By investigating the chemical potential drops in 
the different regions and at the contacts we obtain the 
equilibrium and nonequilibrium Seebeck coefficients, 



_ {Si+ SslPcrl) Li/ai + {Scl + SsclPsi) /Scl + SnLn/cTN + {Sc2 + Ssc2Pt.2) /Sc2 + (52 + Ss2Pa2) L21 02 , , 
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FIG. 15: (Color online) Profiles of the spin potential (a), the 
total chemical potential (b), and the spin current (c) for a 
NigiFeig/Cu/NigiFeig junction at T = 300 K with Li = L2 = 
100 nm, Liv = 50 nm, and AT = -100 mK. The sohd hues 
show the profiles for the parallel configuration, the dashed 
lines for the antiparallel configuration. 



Once more, Eq. (100) has been derived in the limit of 
Li ^ Asi and L2 S> As2, which usually applies to most 
devices. 

Figure 15 shows the profiles for a symmetric F/N/F 
junction consisting of NigiFeig as ferromagnets and Cu 
as the nonmagnetic material (for the corresponding pa- 
rameters see Sees. Ill and IV) for T = (Ti -f-T2)/2 = 300 
K and AT = T2-T1 = -100 mK. Here the lengths of the 
individual constituents are chosen to be Li = L2 = 100 
nm and Xat = 50 nm. The contact parameters are 
= i?^2 = 1 X 10-i<^ f7m2, 5*^ = Sc2 = -1.0 X 10-6 
V/K, and Pei = ±Pe2 = 0.5 and Ssci = ±Ssc2 = O.SS'ci 
depending on whether the parallel (+) or antiparallel (— ) 
configuration is investigated. As shown in Figs. 15 (a) 
and (c), spin is injected into the N region from both F 



regions in the antiparallel configuration. If the F/N/F 
junction is in the parallel configuration spin is injected 
into the N region from one F region, while at the opposite 
interface spin is extracted from the N region. Changing 
the sign of AT would lead to spin extraction from the 
N region in the antiparallel configuration, whereas spin 
would still be injected at one interface and extracted at 
the other interface. In Fig. 15 (b) one can observe a drop 
of the total chemical potential across the F/N/F junction 
for both, the parallel and antiparallel configurations. If 
an asymmetric F/N/F junction (for example by choos- 
ing different lengths Li and L2 or different materials for 
Fi and F2) is considered, the qualitative properties of 
Fig. 15 will remain the same, although the graphs will be 
distorted compared to the symmetric case. 

Next, we look at the difference between the drops of the 
chemical potential [given by Eq. (97)] in the parallel and 
antiparallel configurations (denoted by the superscripts 
i ^ttiti in the following), as a quantitative measure of 
the spin accumulation in the N region (thermal analog 
of the giant magnetoresistance) . If one analyses the tem- 
perature profile T{x) and the local equilibrium chemical 
potential 77 [T(a;)], one finds that within our model they 
are the same for the parallel and antiparallel configura- 
tions (in the limit Sxj <C ^fC). Hence, the difference 
between the drops of the chemical potential is just the 
drop of the quasichcmical potentials, that is. 



A(77/e + /i)" - A(?7/e + ^)i^ = A^f^- A/i^^ (101) 

where A/i* — ^"^{Lfq + L2) — /i*(— Li). Moreover, the 
equilibrium Seebeck coefficients given by Eq. (99) are the 
same for both configurations and consequently 

- Am^^ = {5S^^ - 5S^^) AT, (102) 

which, in the limit of Li ^ A^i and L2 ^ As2, yields 



A/x^^ - A^ 



[(5,iA,i/ai + Sscl/^cl) {R2P2 + i?c2Ps2) + {Ss2Xs2/<y2 + Ssc2/^c2) {RlPl + RclP^l)] A.AtVT, 



N 



Do sinh {Ln/\sn) 



(103) 



if Eq. (100) is inserted for each of the noncquilibrium See- 
beck coefficients. In Eq. (103) as well as in the following 
we choose to express the system parameters in terms of 
the parallel configuration (for example, P2 = Pj^ etc.). 
As mentioned before, in our approximation the temper- 
ature gradient in the N region, given by Eq. (93), does 
not depend on whether the system is in its parallel or 
antiparallel configuration. 

The charge neutrality condition (17) enables us to re- 
late A/Lt' to the voltage drop measured across the junc- 
tion, A(p' = ip^{LN + L2) — ((C*(— Li). Using this the dif- 



ference between the voltage drops in both configurations 
can be written as 



^11 



n _M _ ,,n\ _ 9s2 
91 ^ ^92 
- (A^^^ - A^i^^) , 

(104) 

where the shorthand notations /i*^ = [i\{—L]) and 
= \i\{Lf^ + L2) have been introduced. For Li >• Asi 
and L2 ^ \s2 the contributions to Eq. (104) origi- 
nating from the spin accumulation at a; = —L\ and 
X = Ljv + L2, mIl ~ /^II ^nd /i^Jj -I- /Xj^, are small com- 
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FIG. 16: (Color online) Difference between the chemical po- 
tential drops of the parallel and antiparallel configurations, 
Afi^^ — A/i^'^, as a function of the length of the N region, 
Liv, for a NisiFeig/Cu/NigiFeig junction at T = 300 K with 
Li = L2 = 100 nm and AT = -100 mK. 



pared to A/i^^ — A^u^^^ and consequently 



(105) 



Thus, one can also measure the difference between the 
quasichemical potential drops electrostatically, namely 
as the difference between the voltage drops across the 
F/N/F j miction. 

Figure 16 shows the dependence of A/i^^ — A/i^^ 
on the length of the N region, Ljv, for a symmetric 
NigiFeig/Cu/NigiFeig junction similar to the one con- 
sidered above (apart from L^v, i?ci) and i?c2 the param- 
eters are the same as in Fig. 15) for the contact resis- 
tances i?ci — Rc2 = 1 X 10~^^ f)m^ and Rd = Rc2 — 
1 X 10"^* VlTe? . With increasing length of the N region 
the amplitude of the voltage difference decreases until, 
for very large N regions with Ln ^ \sn, there is no 
difference between the voltage drops in the parallel and 
antiparallel configurations and A/i^^ — A/i^^ — >■ 0. If 
Ljv is comparable or even smaller than the spin diffu- 
sion length {XsN ~ 350 nm in Cu), the voltage drops 
across the F/N/F junction are different for the different 
configurations with A/i^^ — A/x^"^ given by Eq. (103). 



B. Peltier effects in F/N/F junctions 

The section on F/N/F junctions is concluded by a brief 
discussion of Peltier effects in such structures in the limit 
of S\j <C VC and L1/2 3> \si/2- 

Figure 17 (a) summarizes the first system considered: 
A charge current j is driven across an isothermal F/N/F 
junction and there is heating/cooling the interfaces. Sim- 
ilarly to Sec. IV C, the electrical spin injection efficiencies 
at the interfaces, Pji = js(0)/j and Pj2 = js{LN)l j, are 
given by the standard model of electrical spin injection 



FIG. 17: (Color online) A schematic illustration of a F/N/F 
junction in the electrical spin injection setup, where (a) refers 
to an isothermal junction and (b) to the situation where 
jq{x) = 0. The fact that in (b) the temperature at one end of 
the junction is not given as an external boundary condition, 
but has to be calculated from the model is implied by "?" . 



and, as described in detail in Ref. 5, read 

P -pO Rfe RNCOth{LN/XsN)+Rcl 



Ri 



+ PjiRpff 



Dn 



R 



N 



(106) 



Do sinh(LAr/AsAr)' 



where Dq is given by Eq. (92) and k, I — 1,2 and k ^ I. 
The effective resistances of the individual F/N junctions, 
Rpj^, are given by Eq. (91) and their electrical spin in- 
jection efficiencies by 



^jk - 



PskRck + PakRk 
Rk + Rck + Rn 



(107) 



As noted above, the electrical spin injection efficiencies 
of a F/N/F junction, Eq. (106), are composed of the 
electrical spin injection efficiencies of the individual F/N 
junctions in the same way the thermal spin injection ef- 
ficiencies, Eq. (89), are composed of the thermal spin 
injection efficiencies of the individual F/N junctions. 

Consequently, the rates of heat production/dissipation 
at contacts Ci and C2 read 



r™i*=j9(o-)-j,(o+) = r,i + r 



9I' 



(108) 



(109) 



and consist of contributions from the conventional Peltier 
effect, 
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FIG. 18: (Color online) Profiles of the heat current for an 
isothermal NisiFeig/Cu/NigiFeig junction at T = 300 K with 
Li = L2 = 100 nm, Ljv = 50 nm, and j = 10^ A/m^ The 
solid line shows the profile for the parallel configuration, the 
dashed line for the antiparallel configuration. 



as well as contributions from the spin Peltier effect, 

TSaiPjij 



r 



TSs2Pj2j 



q2 



(112) 



(113) 



Figure 18 illustrates this situation for an isothermal 
NigiFeig/Cu/NigiFeig junction (in parallel and antipar- 
allel configurations) at T = 300 K with Li = L2 = 100 
nm, Ln = 50 nm, Rd — Rc2 = 1 x 10~^^ flm^, 
Scl = 5-^2 - -1.0 X 10^6 y/K, = ±p^^ = 0.5, 

Ssci = ±Ssc2 = O.bSci, and j = 10^ A/m^. The pro- 
files of the heat current in Fig. 18 show that — for the 
parameters chosen — there is cooling at Ci {x = 0) as 
heat flows away from it, while heat flows to C2 and leads 



to heating in the region around the C2 {x = Lm)- The 
widths of those regions of heating/cooling are given by 
the individual spin diffusion lengths. 

The second system considered is a F/N/F junction 
where jq{x) = and across which an electric current 
j is driven and one end of which is anchored at a fixed 
temperature [see Fig. 11 (b)]. Requiring the charge, spin, 
and heat currents given by Eqs. (43)-(45) to be contin- 
uous and imposing the additional boundary conditions 
lim fisix) = 0, we find that the temperature drop 

a:— f ±C30 

across the junction, AT = ATch -I- AT^, is composed of a 
drop due to the conventional Peltier effect. 



ATeh = 



{Si + SslPal) Li Scl + SsclPsi SnL 



N 



2C(7i 2/3I]cl IClJjq 

Sc2 + Ssc2Pt,2 {S2 + Ss2Pa2) ^2 



2£f72 



(114) 

and a contribution due to the spin accumulation in the 
region around the interfaces. 



AT, J-^^^^,. (0-) - ^^^i^^M. 



2C 



+ ^^^4^^. (0+)-A^. (0-)] 



2C 

Ssc2 (l — -^22) 

2C 



(115) 

Here we are mainly interested in the difference between 
those temperature drops in configurations of parallel and 
antiparallel magnetizations of the ferromagnets (denoted 
by the superscripts i =tt:t4- as in Sec. VA). With the 
temperature drop due to the conventional Peltier effect 
being the same for both configurations, this difference is 
exclusively due to the spin accumulation, that is, AT^^ — 
AT'^^ = ATP - AT^^, which can be calculated as 



AT'^t _ AT1^^ = KSslXsl/'Jl + Sscl/^cl) {R2P2 + Rc2P^2) + {Ss2K2/<J2 + Ssc2/^c2) {RlPl + RclPsi)] Rwj 

CDosinh{LN/XsN) 

(116) 



r 



where we have expressed the system parameters in terms 
of the parallel configuration (see Sec. VA). 

For illustration the temperature profiles of a 
NisiFeig/Cu/NisiFeig junction at T = 300 K with 
jq{x) — 0, Li — L2 — 100 nm, Ljv = 50 nm, Rd = 
i?c2 = 1 X 10-16 f7ni2, Scl = Sc2 = -1-0 X 10-6 V/K, 
Pei = ±Pe2 = 0.5, Ssci = ±Ssc2 = 0.55,1, and j = 10" 
A/m^ are shown in Fig. 18 (a) for both, parallel and an- 
tiparallel magnetizations in the ferromagnets. While the 
main (linear) contribution to the temperature drop orig- 



inates from the charge Peltier effect and is the same for 
both configurations, the spin accumulation near the in- 
terfaces is different for each configuration and accounts 
for different temperature profiles. Figure 18 (b), which 
depicts the difference between the temperature profiles 
of the parallel and antiparallel configurations, also shows 
that this difference in the temperature profiles arises in 
the F regions near the interfaces and within the spin dif- 
fusion lengths. Outside these regions the temperature 
difference remains constant. 
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FIG. 19: (Color online) Temperature profile (a) of a 
NigiFeig/Cu/NisiFeig junction witfi with j,(a:) = 0, Ti = 300 
K, Li = L2 = 100 nm, Ln = 50 nm, and j = lO" A/m^ The 
solid line shows the profile for the parallel configuration, the 
dashed line for the antiparallel configuration. The profile of 
the temperature difference between the parallel and antipar- 
allel configurations is shown in the inset (b). 



where we found that only at the boundaries of the ferro- 
magnet there is significant electronic spin accumulation, 
which, however, decays within the spin diffusion length 
and can therefore not be responsible for the linear inverse 
spin Hall voltage measured by Uchida ct al.^*" Further- 
more, we have analyzed F/N and F/N/F junctions. For 
F/N junctions we have shown that a temperature differ- 
ence between both ends of the junction generates pure 
spin currents which can be used to extract or inject spin 
at the interface between the F and N regions. We have 
also derived a formula to measure the efficiency of the 
spin injection(extraction). In the case of a F/N/F junc- 
tion a temperature difference can also be used to extract 
or inject spin into the N region if the junction is in a 
antiparallel configuration. Moreover, a formula has been 
derived to calculate the difference between the voltage 
drops across the junction in the parallel and antiparallel 
configurations. Finally, we have investigated the Peltier 
and spin Peltier effects in F/N and F/N/F junctions and 
derived analytical formulas to describe their respective 
contributions to the heating/cooling in these systems. 



VI. CONCLUSION 

We have generalized the standard model of spin injec- 
tion as explained in Refs. 1,2,5 to describe the coupling 
between charge, spin, and heat transport in metals. The 
formalism has then been used to describe the electronic 
contribution to the spin Seebeck effect in such materials. 
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The chemical potential is not only a function of the tem- 
perature but also of the total electron density. 
In general, fj,x also depends on the temperature T. If we 
consider different, space dependent local equilibrium tem- 
peratures T{x), the gradient of the quasichcmical potential 
reads V^a [x, T{x)] = ^ + ^ VT. Since we are only in- 
terested in first order effects, the temperature dependence 
of nx, which leads to a second order contribution (in the 
nonequilibrium quantities nx [x,T{x)], fix), and VT{x)), 
can bo omitted. 

The Einstein relation is obtained by requiring that jx = 
if V I +ij,x{x)^=0 and VT{x) = 0. 

^® In non-degenerate semiconductors one can relax this con- 
dition and obtain the resulting nonlinear current-voltage 
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^"^ This is due to the assumption of constant temperature gra- 
dients and would not be the case if the full model was used. 

^® Here we use that VTn = {aF/crN)'VTF. 



